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O |- Abstract. We discuss the production of hadrons in e + e collisions at y/s = 91 GeV. 

1^ We address the question wether the particle yields measured in the final states are 

' consistent with the statistical model predictions. In the model formulation we account 

for exact conservation of all relevant quantum numbers using the canonical description 
of the partition function. Within our model the validity of the thermodynamical 
approach to quantify particle production in e + e _ annihilations is not obvious. 
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1. Introduction 

CO 
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One of the essential results in heavy ion collisions was the observation that particle yields 
measured in a final state closely resemble a thermal equilibrium population [H EJ [3]. 
The natural question was whether this statistical behavior is a unique feature of high 
energy nucleus-nucleus collisions or whether it is also applicable in elementary collisions 
like, e.g., e + e~. Previous publications [4] indicated that indeed hadron production 
in e + e~ collisions can be well described within a thermal model provided that local 
quantum number conservation is properly implemented. In view of the most complete 
and extended data summarized e.g. by the Particle Data Group (PDG) [5] the above 
question has been recently addressed independently in [6] and [7]. 

In this contribution we discuss hadron production in e + e~ annihilations at ^/s = 
91 GeV based on the thermal model analysis of [6j. We also focus on different 
implementations of the statistical model and discuss the importance of quantum statistic 
effects and the mass-cut in the hadron mass spectrum. We quantify the production of 
heavy flavors and compare the model predictions with available data. 
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2. The statistical model and charge conservation 



The usual form of the statistical model in the grand canonical ensemble formalism cannot 
be used when the number of produced charged particles is small. This is the case if either 
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Figure 1. The (left-hand figure) shows deviations of pion and kaon yields from 
their exact quantum statistics values (see text). The k = 1 term corresponds to the 
Boltzmann approximation. The (right-hand figure): the relative change of the particle 
yields calculated in the hadron resonance gas model with a mass spectrum cut at the 
mass M = 1.7 and M = 3.0 GeV. 



the temperature T or the volume V or both are small. As a rule of thumb one needs 
VT 3 >> 1 for a grand canonical description to be applicable [$]. In e + e~ annihilations 
where modelling the particle production within a thermal approach requires an exact 
formulation of conservation laws. In such a system one needs to account for an exact 
conservation of five quantum numbers: the baryon number N, strangeness S, electric 
charge Q, charm C and bottom B. 

The appropriate tool to deal in a statistical mechanics framework with a system of 
quantum numbers X = (N, S, Q, C, B) is the canonical partition function [HI [9] 

Zn,s,q,c,b = j^y 5 j d ^ ^ m ex P (E *i)> W 

where 

Zj = g 3 j^ J d 3 p ln(l ± exp (-^ + m 2 /T _ ^.0))±i, (2) 

and Xj is a five component vector Xj = (Nj,Sj,Qj,Cj,Bj) containing the quantum 
numbers of the particle species j. The quantity (f> = (0jv> 4>s, 4>Qi (fic, (Pb) is an element of 
the symmetry group [[/(l)] 5 related with additive conservation laws. In this expression, 
each (px corresponds to the conservation of the corresponding quantum number X and Zj 
is the single particle partition function for particle with mass mj, spin-isospin degeneracy 
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factor gj, and a system with volume V and temperature T. The sum in Eq. ([T]) runs 
over all particle species in the hadronic gas. 

The integral representation of the partition function in Eq. (TjQ) is not convenient 
for numerical analysis as the integrand is a strongly oscillating function. Thus, we first 
expand the logarithm 

OO 

^±^ = £(±1)^ (3) 

k=l 

and then, using the method described in [TO], [Til 12], we express the partition function 
(PQ) in a series of Bessel functions to obtain a result that is free from oscillations. 
Furthermore, from Eq. (JTJ) we obtain the multiplicity (rij) for particle species j by 
introducing a fugacity parameter Xj which multiplies the particle partition function 
Zj and by differentiating 

d\nZ 



( nj ) 



d\ 3 



(4) 

A,=l 



The first term in the expansion of the logarithm in Eq. ([3]) corresponds to the 
Boltzmann approximation which is suited only for nij » T. Such a condition is 
satisfied for baryons since their masses are larger than a typical temperature of the 
hadron resonance gas which never exceeds a critical value T c ~ 200 MeV required for 
deconfinement. However, for light bosons like pions or kaons the quantum statistics 
is of importance as the temperature is comparable to their masses. Fig. (JTjTeft) 
shows relative deviations of pion and kaon multiplicities from their quantum statistics 
values with increasing numbers of terms k in the expansion ([3]). The calculations were 
performed for T = 157 and V = 32 fm 3 . It is clear that the Boltzmann approximation 
is by far not sufficient to reproduce the quantum statistics results. The pion yield under 
Boltzmann approximation deviates by more than 7% from the exact quantum statistics 
result. For kaons this difference is only 1% due to the larger mass. For pions, five 
terms are needed in the expansion (J3]) to get quantum statistics value with deviations 
below 10~ 3 . For heavier particles, for instance for protons, the quantum and Boltzmann 
statistics differs by less than 0.1%. Thus, in the model comparison with e + e~ data 
we apply quantum statistics for light mesons and we use Boltzmann statistics for all 
baryons. 



3. Modelling the e + e events within the statistical approach 

The hadron multiplicity calculation within the statistical model basically proceeds in 
two steps. First, a primary hadron yield Nff 1 , is calculated using (CQ) and (0J. A crucial 
assumption of the model is that the final yields of all particles are fixed at a common 
temperature, the chemical decoupling point. As a second step, all resonances in the 
gas which are unstable against strong decays are allowed to decay into lighter stable 
hadrons, using appropriate branching ratios (B) and multiplicities (M) for the decay 
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Figure 2. Comparison between the best fit thermal model calculations and 
experimental hadron multiplicities (sum of particles and antiparticle yields and for 
the 2 jets) for e + e~ collisions at y/s=91 GeV. The upper panel shows the fit in an 
uncorrelated jet scheme when including the feed-down contribution from heavy quarks. 
The middle panel is after subtraction of this contribution. The results in the lower 
panel are obtained in a correlated (5-flevour) jet scheme fit where the flavor abundances 
are extra input parameters from data (see text). The best thermal model fit parameters 
are listed for each case. 



j — > h published by the PDG [13]. The abundances in the final state are thus determined 
by 

N h = Nl h + J2 ■ B(j - h)M(j - h) (5) 

3 

where the sum runs over all resonance species. 

From Eq. (jSJ) it is clear, that the final multiplicity of stable hadrons depends on 
the number of resonances used in the sum. In general one should include contributions 
of all known resonances as listed by the PDG. Fig. (Tjjright) shows the relative change 
of different particle yields when applying the mass cut M = 1.7 and M = 3.0 GeV in 
Eq. (jSJ). For mesons this difference amounts to 5% and is as large as 15% for baryons. 
Thus, it is clear that restricting the mass spectrum only to the resonances with the 
mass M < 1.7 GeV might be not sufficient at the level af accuracy of data of a few 
percent as is the case in e + e~ collisions. However, our knowledge of decay properties of 
heavier resonances is by far not complete, which causes systematic uncertainties of the 
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statistical model. The importance of heavier resonances originating from the Hagedorn 
mass spectrum in the analysis of particle production in heavy ion collisions has been 
recently analyzed in [H] . In the following we account for resonance contributions up to 
the mass M < 3 GeV. 

In general, the resonance gas model formulated in the canonical ensemble is 
described by only two basic thermal parameters: the temperature T and the volume 
V of the system. To explore a possible strangeness undersaturation we introduce an 
additional parameter 7 S into the partition function to account for a possible deviation 
of strange particle yields from their chemical equilibrium values [I]. When applying 
the statistical model to particle production in e + e~ annihilations we have to take into 
account that most hadronic events in such collisions are two-jet events, originating 
from quark-antiquark pairs of the five lightest flavors. Since we would like to address 
the issue of overall equilibration in these systems, we have to specify how the initial 
quantum numbers are distributed between the two-jets. We will consider two scenarios: 
an uncorrelated and correlated jet scheme. 

In an uncorrelated jet scheme each jet is treated as a fireball with vanishing quantum 
numbers as fixed by the entrance channel. It is clear at this point that hadrons from 
jets with heavy quarks (c and b) will be greatly underestimated by the model because 
of the large Boltzmann suppression factors. In this approach the issue of equilibration 
is effectively addressed only for hadrons with light quarks (u, d, s). It is important to 
recognize that the measured yields of these hadrons contain the contribution from the 
e + e~ annihilation events into cc and bo. Heavy-quark production is indeed significant 
and is very precisely measured, in particular at the Z mass, where the measurements 
are very well described by the standard model . Hence, heavy-quark production is 
manifestly non-thermal in origin. We therefore consider two cases: i) we fit the data 
as measured and ii) we subtract from the yields of hadrons carrying light quarks the 
contribution originating from charm and bottom decays based on available data for the 
charmed and bottom hadron production and their branching ratios. 

In a correlated jet scheme (case (iii) in Fig. ([2]) ) the initial quantum numbers are 
distributed such that each jet carries quantum numbers of either S — ±1, C — ±1, 
B = ±1, or vanishing quantum numbers in the case of uu and dd jets. The fractions 
of the quark flavors in hadronic events [5] are external input values, unrelated to the 
thermal model (see also Table II in ref. [7]). 

3.1. Model comparison with e + e~ data at ^/s = 91 GeV 

For the fit procedure we use the complete set of all measured hadron yields with 
the exception of those containing charm or bottom quarks. A x 2 fit is performed by 
minimizing 




(6) 
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as a function of the three parameters T, V and 7 S , taking account of the experimental 
uncertainties cr h . 

The resulting best fit to the data for a correlated and uncorrelated jet scheme 
is shown in Fig. (J2J). The two cases of the fit, without and with the subtraction of 
the contribution from heavy quarks in an uncorrelated jet scheme are also shown in 
this figure. We first note the overall behavior of the data, namely an approximately 
exponential decrease of the particle yield with an increasing particle mass. Such a 
behavior is expected in the hadron resonance gas model due to the Boltzmann factors, 
thus indicating the presence of statistical features of hadron production in elementary 
collisions. The quantitative description of the data with the statistical model is, however, 
rather poor and certainly no improvement is visible for the case of subtracting charm 
and bottom contributions. The low fit quality is reflected through the large x 2 values 
per degree of freedom. In addition, discrepancies between individual data points and 
fit values larger than 5 standard deviations are not rare. There are also problems in 
determination of the fit parameters: T, V and 7 S . In \ 2 contour plots, both in (T, V)- 
and (T, 7 s )-plane, one notices strong anticorrelations between fit parameters [6]. In 
addition there is a series of local minima which makes it difficult to uniquely determine 
the model parameters. Such local minima are typical for poor fits and imply large 
uncertainties in the determinations of the fit parameters. 

3.2. Model description of heavy quark hadron production 

In the canonical formulation of the statistical model with exact charge^! conservation 
the abundances of charged particles depend crucially on the overall charge in a system. 
In order to illustrate this let us consider a model where only one charge, e.g. charm, 
is conserved exactly. In such a case the multiplicity < N,- t >c-=±i °f particle i with 
mass rrii that carries charge d = ±1 in a system of the total charge C, volume V and 
temperature T is obtained under the Boltzmann approximation from [8]: 

<N . >C y Cl Z Tl I CT i(2Vx) 

<^> Ci=±1 -V Zl x (7) 



where 



z, 



' --^m*TK 2 (m t /T), Z ±1 = Vzf =± \ x = y/Z^l, (8) 



and where Ik and K 2 are Bessel functions and the argument of quantifies the total 
number of charged particle pairs. 

From Eq. ((Zj) one recognizes an essential difference in particle yields if the total 
charge C = or C = ±1. Indeed, assuming that Vx < 1 and applying an asymptotic 
expansion of the Bessel functions one finds e.g. that multiplicities of particles with 
charge Cj = +1 is 

< Ni >~ V 2 Zi Z_ 1: (9) 

| The charge is consider here to be any quantum number related with U(l) symmetry. 
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if the total charge of a system C = 0, and 

z +i 

if the total charge of a system C — +1. 



(10) 
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Figure 3. Comparison between the thermal model calculations in a correlated jet 
scheme and experimental data on charmed and bottomed hadrons for e + e~ collisions 
at ^=91 GeV. Also shown are the model parameters used in the calculations. 



For a charge neutral system the charged particles yield ([HD are strongly suppressed 
due to canonical effects, since a particle has to be produced in a pair with an antiparticle 
in order to fulfill charge neutrality. This is a well known "canonical suppression effect" 
which is crucial e.g. to quantify "strangeness enhancement" and production in hadron- 
hadron and heavy ion collisions at lower energies [2]. For the total charge C = ±1 
the canonical effect results in enhancement of particle yields as shown in Eq. ([TO]) . 
The charge C — ±1 of a system is here redistributed between all particles that carry 
Ci = ±1 with weights given by the ratio of the thermal phase-space of particle % to the 
phase-space of all negatively or positively charged particles. Here, even the Boltzmann 
suppression is to a large extent cancelled out in the ratio (see Eq. ( 1101) ). 

The above examples imply a qualitative difference in open charm and bottom 
production in e + e~ collisions when using an uncorrelated and correlated jet scheme 
scenario. If C = B = in a jet then the thermal phase-space of charm and bottom is 
strongly suppressed as in Eq. ([9]), leading to large discrepancies between model results 
and data. One needs to check whether in the correlated jet scheme, where the open 
charm and bottom are statistically enhanced as in Eq. (Q, the measured yields are 
comparable with model predictions. 

In quantitative analysis, one calculates the heavy quark yields with the partition 
function (JTj) rather than with the approximate Eq. (ITUl) . Fig. (j3J) shows the statistical 
model results for charmed and bottom particles obtained in the correlated jet scheme. 
One recognizes a good description of data by the canonical statistical model. Thus, if the 
fireballs related with each jet are charged then the model description of the distribution 
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of these charges between different particles agrees with data. The value of x 2 l^df — 2 
indicates a good quality of the model description of data. We have to stress that, in the 
canonical model, yields of hidden bottom and charm meson are still strongly suppressed 
by the Boltzmann factors, thus their multiplicities deviate from data. Indeed, in Fig. 
(J3J), the yield of Y is larger by several orders of magnitudes than the model results (not 
shown in Fig.(j3])), indicating a non-thermal origin of this particle. This is not the case 
for ip since in e + e~ annihilations at y/s = 91 GeV as they are almost entirely originating 
from decays of bottom [T5]J§I. 



4. Summary and conclusions 



We have analyzed the experimental data on hadron yields with light and heavy quarks 
in e + e~ collisions at \fs = 91 GeV within the statistical model. The conservation of five 
quantum numbers was included in the framework of the canonical partition function. 
Our results qualitatively confirmed the previous finding [I] that statistical features 
are present in hadron production in e + e~ annihilation. The resulting temperature of 
160 — 170 MeV lies in the bulk expected at the chemical freezeout in heavy ion collisions 
at high energy [U [2] and agrees with first results obtained in e + e~ systems Our 
results on open charm and bottom hadrons, using the measured c and b fractions of 
jets, showing good agreement with data also confirmed earlier observation [I] . However, 
in view of a rather poor fit to measured yields of hadrons with light quarks and a clearly 
non-thermal origin of the hidden charm and bottom particles as well as essentially 
different characteristics of the collision fireball in e + e~ and in heavy ion collisions [6], 
the general validity of the thermodynamical approach to the particle production in e + e~ 
annihilation at LEP energies is not obvious. 
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